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The properties of M06S3I6 nanowires were investigated with ab initio calculations based on the 
density-functional theory. The molecules build weakly coupled one-dimensional chains with three 
sulfur atoms in the bridging planes between the Mo octahedra, each dressed with six iodines. Upon 
uniaxial strain along the wires, each bridging plane shows two energy minima, one in the ground 
state with the calculated Young modulus Y — 82 GPa, and one in the stretched state with Y — 94 
GPa. Both values are at least four times smaller than the experimental values and the origin of the 
discrepancy remains a puzzle. The ideal tensile strength is about 8.4 GPa, the chains break in the 
Mo-Mo bonds within the octahedra and not in the S bridges. The charge-carrier conductivity is 
strongly anisotropic and the M06S3I6 nanowires behave as quasi-one-dimensional conductors in the 
whole range of investigated strains. The conductivity is extremely sensitive to strain, making this 
material very suitable for stain gauges. Very clean nanowires with good contacts may be expected 
to behave as ballistic quantum wires over lengths of several fim. On the other hand, with high- 
impedance contacts they are good candidates for the observation of Luttinger liquid behaviour. The 
pronounced ID nature of the M06S3I6 nanowires makes them a uniquely versatile and user- friendly 
system for the investigation of ID physics. 

PACS numbers: 61.46.-w, 62.25.+g, 73.22.-f, 71. 15. Mb 



I. INTRODUCTION 

Inorganic nanowires and nanotubes are rapidly gaining 
in importance as materials of great scientific and techno- 
logical interest because they exhibit a substantially en- 
hanced level of functionality which can open novel oppor- 
tunities for devices and applications. The main advan- 
tage of inorganic nanowires over other related materials 
is the possibility of controlling the relevant physical prop- 
erties by selective engineering of their geometry and/or 
composition. Molybdenum-based inorganic nanowires 
and nanotubes in particular are emerging as materials 
with very promising physical properties because they are 
easier to functionalize and to synthesize in clean, single 
stoichiometry form than other one-dimensional materi- 
als. The synthesis, some physical properties and appli- 
cations of several molybdenum-based nanowires have al- 
ready been reportediSii^ Recently, significant research 
effort has been focused on nanowires composed of Mo, 
chalcogens (S) and halogens (I) in the form Mc>6Sg_ x I x 
(MoSIa;)A£ which are best described as molybdenum 
chalcogenide-halide polymers. They have a backbone 
skeleton composed of M06 octahedra, each dressed with 
six anions (S or I), bound together into one-dimensional 
chains by three bridging anions (either S or I). A spe- 
cific feature of these materials is the growth of identical 
chains in the form of bundles. The materials have strong 
anisotropy, large Young moduli along the wires, very 
small shear moduli, and controllable electronic proper- 
ties. In contrast to related nanowires such as I^Mo^Sc^ 
which decompose rapidly in air, the MoSIx nanowires are 
entirely air-stable. Yang et ali have recently investigated 
the M012S9I9 nanowires theoretically and have pointed 
out the bistable nature of the sulfur bridges, which allow 
the wires to stretch in an accordion-like fashion at vir- 



tually no energy cost. They also pointed that strain can 
have a significant effect on the electronic band structure, 
and, interestingly, the appearance of a net spin polarisa- 
tion within the nanowires in the unstrained position. 



In this paper we report on the theoretical investiga- 
tions of extraordinary elastic and electronic properties 
of the most common reproducibly synthesized MoSI^ 
nanowire material, namely M06S3I6. Meden et al£ re- 
ported the structure where, after purification, the X-ray 
diffraction experiments indicated an iodine occupancy 
of the bridging atoms and only a small degree (~ 10 
%) of disorder in the anion occupancies (either S or I). 
More recently, Nicolosi et al£ reported on the electron 
microscopy of same isomer with the same skeletal struc- 
ture but with the bridging planes occupied by the sulfur 
atoms and all the iodines attached to the Mo octahedra. 
As we will comment later, both isomers have almost the 
same energy therefore it is possible that both are syn- 
thesized at similar rates. The question arises, how to 
identify them, what are the similarities and what are the 
differences between the two isomers. We shall concen- 
trate in this paper on the isomer with three S atoms in 
the bridging planes, Fig. ^ From the tribology experi- 
ments we know that the interaction between individual 
nanowires is weak, 9 therefore we focus our analysis on 
the longitudinal properties of individual nanowires, or- 
dered in a hexagonal array according to the spacegroup 
R3c. This space group assumes the inversion symmetry 
which makes all the sulfur bridges equivalent. As a con- 
sequence we can scan the whole range of strains. Using 
the density-functional theory (DFT) we investigate the 
effect of strain on the elastic and electronic properties of 
these nanowires. 
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FIG. 1: (Color online) Top (upper panel) and side (bottom 
panel) view on the primitive cell of M06S3I6 nanowires or- 
dered according to the spacegroup R3c. Mo octahedra (grey) 
are dressed with 6 I atoms (violet - dark grey) and separated 
by 3 bridging S atoms (yellow - light grey) . The X-ray scatter- 
ing data and the DFT calculations of the R3 isomer indicate 
a ~ 6° rotation of the chains around the hexagonal c axis and 
a certain degree of occupational (either S or I) disorder. The 
distant I atoms belong to the neighbouring chains. 



II. COMPUTATIONAL DETAILS 



properties of strained crystals. We thus obtain the elas- 
tic constants en which are close to the Young moduli Y 
because of weak interchain coupling. 



III. RESULTS 

The DFT simulations of Mo-S-I arranged as shown in 
Fig. n gi ye the equilibrium hexagonal lattice constants 
c = 11.56 A and a = 16.54 A which agree to within 3% 
with the experimental values c = 11.95 A and a = 16.39 
A although we assumed a different spacegroup (the simu- 
lated group was R3c, while the X-ray diffraction suggests 
P63 or lower) & In this configuration the Mo-S-Mo angle 
across the bridging plane is 90.8°. Upon stretching the 
nanowires we find another energy minimum at c = 13.76 
A with the Mo-S-Mo angles 145°, Fig.^a). Yang et aU 
ascribed the two equilibrium configurations to different 
s-p hybridizations of sulfur. However, inspection of the 
atomically resolved electron density of states, discussed 
later, shows very little s-p hybridization on sulfur. On 
the other hand, in the second energy minimum, the sep- 
aration between the S atoms in the same bridging plane, 
ds-s — 3.79 A, is close to twice the van der Waals radius 
(r v dw = 1-85 A) or to the separation between S in neigh- 
bouring layers in layered M0S2, d vc iw = 3.49 A. There- 
fore we suggest the possibility that the van der Waals 
interaction between S on the same bridging plane stabi- 
lizes the second energy minimum. 

The deeper minimum (unstrained wire) is subject to a 
potential barrier AE = 0.47 eV and the higher one by 
AE = 0.26 eV, both per formula unit. These barriers 
are too high to be thermally excited, the only way to 
move from one to the other minimum, is by applying an 
external stress. 



The structure and properties of Mo-S-I nanowires 
were analysed using the WIEN2k code. 1 — The simu- 
lations were performed on a mixed basis set of aug- 
mented plane waves plus local orbitals (APW+lo)Ai for 
low orbital momenta (t < 2) and linearized augmented 
plane waves (LAPW) for all the higher orbital momenta. 
The exchange and correlation potential was treated in 
the Perdew, Burke and Ernzerhof generalized-gradient 
approximation^ The muffin-tin radii were set to -R^o = 
1.06 A for Mo, Rf T = 1.27 A for S, and i?j MT = 1.38 A 
for I, the kinetic energy cutoff was E^ ax = 12.3 Ry and 
the plane- wave expansion cutoff was E^ x = 196 Ry. 

First the hexagonal lattice constants a = b, c and 
the atomic coordinates were optimized. In this stage 
the energy was calculated on a tetrahedral mesh with 
44 Appoints in the irreducible part of the Brillouin zone. 
The same density of fc-points was used for the analysis 
of elastic properties. Later, when calculating the elec- 
tronic density of states, 189 fc-points were used. Since 
the interactions between neighbouring chains are weak, 
we kept the lattice constants a and b fixed, varied c and 
optimised the atomic coordinates when calculating the 



A. Elastic properties 

Like other MoSI^ nanowires, also the M06S3I6 
nanowires have large interchain separations with weak, 
van der Waals coupling between the chains. As a conse- 
quence, the nanowires are very anisotropic in their prop- 
erties, including elasticity. They have large elastic con- 
stants in the direction of the wires, en, and small shear 
and Young moduli perpendicular to the wires. Here we 
shall concentrate on the most interesting elastic constant, 

We investigated the deformation of the M06S3I6 chains 
upon longitudinal tensile stress by keeping the lateral lat- 
tice constant o constant, to see how do the chains stretch 
and where are the weakest points in the chains. The 
corresponding calculated stress-strain diagram is shown 
in Fig. |2Jb). Initially, most of the strain is connected 
to the deformation of the Mo-S-Mo bond angles across 
the bridging plane like in an accordion while the Mo-S 
distances and the structure of the Mo octahedra remain 
almost constant, (see Fig. [3J. In this region the Young 
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FIG. 2: (Color online) (a) Strain dependence of the binding 
energy per primitive cell (two M06S3I6 molecules). The two 
minima correspond to short and long Mo-S-Mo bridge con- 
figurations, respectively. In the unstable region the wire has 
short and long Mo-S-Mo bridges simultaneously. The inserts 
show the Mo-S-Mo bridges in the two energy minima. The 
straight line indicates the critical stress, i.e., the equilibrium 
transition between the two minima, (b) Stress-strain graph of 
MoeSale nanowires. The slopes of straight lines indicate the 
elastic constants of the two stable configurations. The max- 
imum at en = 0.315 is the ideal tensile strength of M06S3I6 
nanowires with all the Mo-S-Mo bridges in the long configu- 
ration. 

modulus is en = 82 GPa. Very soon the value of the 
critical stress, a c = 0.39 GPa is reached, (see Fig. Ufa)), 
when it becomes energetically more favourable for the 
Mo-S-Mo bonds to jump into their stretched configura- 
tions. However, the potential barrier between the two 
configurations prevents such transitions, the system re- 
mains in a metastable state (dashed line in Fig. [2J un- 
til the unstable region is reached when all the Mo-S-Mo 
bonds suddenly jump at constant stress to the stretched 
state. On contraction, upon reducing the stress, the Mo- 
S-Mo bonds remain in the stretched state down into the 
metastable region until they eventually jump into the 
slightly compressed state in the unstretched state. The 
M06S3I6 nanowires thus show a hysteresis in the stress- 
strain graph, caused by the potential barrier between the 
two energy minima which originate in the nature of the 
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FIG. 3: (Color online) Variation of interatomic distances d 
and Mo-S-Mo angle a upon stretching of MogSsIs nanowires. 
Vertical dashed line shows the position of the second energy 
minimum. Mo-Mo (o) and (b) denote the Mo-Mo distances 
in the octahedron and across the bridge, respectively. The 
weakest segments where the wires break are the Mo octahe- 
dra. 



Mo-S-Mo bonds across the bridging planes. 

Of course, the Mo-S-Mo angles cannot open beyond 
180°, and this is the reason for an increase in an at 
large strains and for a higher elastic constant, en =94 
GPa, than in the first energy minimum. The non-linear 
regime at higher en is the consequence of anharmonic 
forces and is not related to plastic deformation which is 
absent in our DFT calculations. Of course, the chains 
eventually break at the maximal stress. This defines the 
ideal tensile strength, cr max « 8.4 GPa. Inspection of 
Fig. El tells us that the weakest links, where the chains 
break, are - surprisingly - the Mo octahedra and not the 
Mo-S-Mo bridges. This is also the reason why all MoSI^ 
nanowires have similar tensile strengths. 

In the experiments the wires were subject to an ini- 
tial stress which most probably brings all the Mo-S-Mo 
bridges into their stretched configurations. Only after 
this initial treatment the elastic constants were mea- 
sured. Therefore one must compare the experimental 
elastic constant with that of stretched Mo-S-Mo bridge 
configurations. However, it remains a mystery, why the 
experimental longitudinal Young moduli exceed the clas- 
tic constants en reported in this paper by a factor of 4 
or more^ The above values of en were calcualted on 
a hexagonal array of molecular chanin with all chains 
carrying equal load. For comparison, the Young modu- 
lus of a hexagonal array of single-wall carbon nanotubes 
(SWCNT) with the lattice spacing 1.7 nm is - 600 GPa 
which is about 6 times the calculated value of 

The accordion effect is present also in other one- 
dimensional carbon chains like polyethylene or poly- 
acetylene. What is unique to M06S4I6 is the double 
energy minimum. An additional sulfur in the centre of 
the bridging planes in M06S4I6 significantly increases the 
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FIG. 5: (Color online) Density of states in the conduction 
(valence) band. For Mo, S and I the partial charges in the 
muffin-tin spheres are plotted, the rest are in the interstitial 
region. The main contributions to the conduction band are 
from Mo-4d, S-3p and I-5p orbitals. The letters at the top 
denote the energy ranges, discussed in the text. 



FIG. 4: Bandstructure of the conduction (valence) band along 
the special lines of the hexagonal Brillouin zone. Three sub- 
bands (one of them doubly degenerate) cross the Fermi energy 
and contribute to the charge carrier transport. 

elastic constant (en = 114 GPa), because it hinders the 
"accordion effect" of the Mo-S-Mo bonds. In case of 
M06S4I6 there is also no sign of a second energy mini- 
mum in strained nanowires. Finally we notice that the 
calculated elastic constant c\\ of M06S3I6 is more than 3 
times smaller than en of MogSg or K 2 Mo6S6 where the 
flexible Mo-S-Mo bridges are absent^ 



B. Electronic properties 

The electron bandstructure of the conduction or va- 
lence band, Fig. shows several interpenetrating sub- 
bands, three of them crossing the Fermi energy. Because 
of the anisotropy in the crystal structure the dispersion 
in the lateral directions is very small, we can classify 
Mo 6 S3l6 as a quasi-one-dimensional conductor. The con- 
duction band consists of hybridized Mo-4g?, S-3p and I-5p 
orbitals. The character of the sub-bands is most clearly 
seen in the electron density of states (DOS), Fig.JSJ Close 
to Ep (region A in Fig. 0, the hybridization of the Mo- 
Ad and S-3p states is the origin of good conduction along 
the Mo-S-I chains. This is not surprising since a conduc- 
tion electron on its way along the chain must pass the 
sulfur bridges (S-3p states). This is also seen in the real- 
space charge-density plot in Fig. E{a) where the charge 
density in the energy interval A is seen to run continu- 



ously throughout the chain. A cross-cut through the Mo 
(0001) plane, Fig.^Jb), reveals a combination of metallic 
bonds (with a uniform distribution of electrons within the 
Mo octahedron) and covalent bonds between neighbour- 
ing Mo atoms (a covalent bond is directional between the 
atoms). The I atoms, which are close to this plane, are 
only very weakly bound to the Mo octahedron in this 
energy window. Lower in energy, i.e., in energy window 
B, a strong Mo-4d double peak, also responsible for the 
Mo-Mo bonding, dominates the DOS. 

Towards the bottom of the conduction band (energy 
interval D) the DOS is dominated by the states belonging 
to the dressing atoms (I-p in our case and S-p in case of 
MogSg 15 ) and to polar covalent bonding between Mo and 
I. 

The regions C and D of the conduction-valence band 
are mainly responsible for the Mo-S and Mo-I bonding. 
We identify these two bonds as polar covalent because of 
partial ionization of I and S. In fact, the charges of the 
muffin-tin spheres are as follows: Each Mo gives away 
3.5 electrons, S 0.9 electrons and I 2 electrons, all these 
electrons fill the interstitial region. For comparison, Mo 
in metallic state gives away (into the conduction band) 
only two electrons. These charges are also in accordance 
with electronegativity. For a Mo-S bond the difference 
in electronegativity is 0.42 on the Pauling scale whereas 
for a Mo-I bond it is 0.50, I is slightly more negatively 
charged than S. 

Comparison with the DOS of MogSg chains which have 
no S bridges^ shows similar occupations of Mo-4g? states 
in both cases but a much weaker occupation of the S-3p 
states in case of MogSe close to Ep. 
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FIG. 6: (Color online) Valence charge density in the energy 
window A on the (1100) plane through the centre of a M06S3I6 
nanowire (a) and on the (0001) plane through Mo atoms 
(b). The isolines and the color coding are in units of 0.0167 
electrons/ A 3 . Clearly seen are the charge concentration on 
Mo octahedra and through the bridging sulfur atoms. 



The bands below the conduction band are localized 
to the respective atoms and have well-defined characters 
with very little hybridization. In the whole DOS we see 
very little s-p hybridization in the sulfur orbitals from 
which we conclude that the s-p hybridization cannot be 
responsible for the stabilization of the second energy min- 
imum at c = 13.7 A in M06S3I.6. 

From Fig. 0] the estimated Fermi velocities are in the 
range between 0.7 and 1.3 xlO 6 m/s. With a scattering 
lifetime r - 0.7 x 10~ 14 s (r = 0.1 eV) this leads to a 
mean free path along the wires of the order 5 to 9 nm, 
enough to justify a delocalized description of the charge 
carriers. 

Comparison with the bandstructure of the MogSsIg 
isomer with iodines in the bridging planes (spacegroup 
R3)£ shows that the bandstructure is very sensitive to 
the anion occupancies, in particular to the occupancy 



of the bridging atomic positions. Although the present 
R3c structure has virtually no calculated energy gap, the 
spacegroup R3 isomer has a ~ 0.6 eV bandgap above 
Ep, According to our DFT calculations both structures 
have very similar formation energies, therefore we expect 
that both are synthesized simultaneously in comparable 
quantities. For the same reason we cannot exclude the 
presence of other MogSsIg nanochains with more disorder 
in the anion occupancies. The bandstructure, the corre- 
sponding optical spectra and the longitudinal conductiv- 
ity can be an effective criterion to distinguish between 
different isomers. 

To gain insight on the conductivity behaviour we cal- 
culate the dielectric tensor e(w) in the random-phase ap- 
proximation and in the limit q — > Oil The static metal- 
lic conductivity is associated with the Drude peak which 
comes from transitions between states close to Ef- From 
the Drude peak damped with r = 0.1 eVwe find the lon- 
gitudinal conductivity mi ~5x 10 3 S/cm. This is a fac- 
tor ~ 4 lower than mi of MogSg nanowires, which we at- 
tribute to the weaker conduction of S bridges. The lateral 
conductivity is one order of magnitude smaller as well. 
This leads to a very short carrier mean free path in the 
lateral directions. Clearly the band picture breaks down 
in the lateral directions and the charge-carrier transport 
is via hopping between neighbouring chains. Of course, in 
real systems with finite-length chains the charge carriers 
must hop also between neighbouring chains and this can 
also hinder the charge carrier transport along the chains. 
In the calculation, F = 0.1 eV was taken arbitrarily - for 
very clean nanowires, F could be substantially less and 
the conductivity correspondingly higher. With F = 0.1 
the longitudinal conductivity of M06S3I6 is about two or- 
der of magnitudes smaller than a reported for SWCNT 
at 300 K which are considered to be one of the best one- 
dimensional conductors. 16 



C. Electronic properties of stretched nanowires 

The bandstructure of MoSI^ nanowires and conse- 
quently also the Fermi surfaces and electrical conductiv- 
ity are extremely sensitive to strain along the hexagonal c 
axis. In Fig.[7|we show the bandstructure along the [001] 
symmetry direction for four different strains. Although 
it is difficult to follow a particular sub-band through the 
whole range of strains, we observe that in general the 
bandwidth of the sub-bands that cross Ep and which are 
responsible for the charge-carrier transport first increases 
with strain up to the second energy minimum (c = 13.7 
A) and then decreases. Not surprisingly, the bandwidth 
is related to the Mo-S distance, Fig. [31 which reaches 
its minimum around c = 13.7 A. Beyond c = 13.7 A 
the Mo-S distance increases and the sub-band width de- 
creases, as expected. Similar behaviour is expected for 
the longitudinal DC conductivity. Above 15.2 A a drops 
rapidly because the chains start to break apart. Assum- 
ing the same F in the lateral direction we find that the 
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c — Co = 0.114 nm c — Co = 0.167 nm 




c — Co = 0.220 nm c — Co = 0.325 nm 



FIG. 7: Bandstructures of M06S3I6 nanowires subjected to 
longitudinal strain. Several sub-bands cross the Fermi energy 
and the width of these sub-bands is maximal close to the 
second energy minimum at c — Co = 0.220 nm. 



anisotropy in the static conductivity is large, typically of 
the order 100 for stretched nanowires. The Fermi surfaces 
of MogSalg in the second energy minimum (c = 13.76 A, 
Fig. [SJ show typical features for a quasi-one-dimensional 
material, almost planar Fermi surfaces with very little 
curvature. The number of Fermi surfaces and the Fermi 
velocities are very sensitive to longitudinal strain, and so 
is the conductivity. 

The M06S3I6 nanowires thus behave as quasi-one- 
dimensional metals in the whole range of strains. This 
is different to the M012S9I9 nanowires where a transi- 
tion from metallic to narrow-gap semiconductor has been 
reportedi* 

The lateral conductivity is dominated by hopping to 
neighbournig chains through the bridging sulfur atoms, 
see FigEl Upon longitudinal strain, the sulfur atoms 
move towards the centre of the wires whereas the in- 
terchain separation, dominated by the non-conducting I 
atoms, remains almost constant. As a consequence, the 
lateral conductivity drops rapidly with tensile strain. 

The variation of the conductivity with strain is de- 
scribed with the strain gauge factor, G = Ap/(pAe) 
where p = 1/a. For the longitudinal conductivity we 
find G of the order -10 which is somewhere in between the 




FIG. 8: Fermi surfaces and the Brillouin zone of the rhombo- 
hedral unit cell of stretched nanowires (c= 13.76 A). 

typical values for semiconductors and for metals. How- 
ever, the unusual feature of MoSI-r nanowires, and great 
potential advantage in applications is their extreme flex- 
ibility (large Ae) arising from the accordion effect in the 
Mo-S-Mo bridges. Values of Aer/er ~ 1 can be reached in 
the elastic regime which is an order of magnitude better 
than Act/ct in classical metallic strain gauges. 

IV. DISCUSSION AND CONCLUSION 

The M06S3I6 nanowires investigated in this paper have 
some similarities to the skeletally equivalent M012S9I9, 
but also some important differences. Both nanowires 
have three sulfur atoms in the bridging planes between 
the Mo octahedra and each bridging plane shows two 
energy minima upon uniaxial strain, one in the ground 
state (unstretched wire) and one corresponding to the 
stretched wire. The potential barrier between the two 
minima is high and prevents thermally excited transitions 
between the two minima, it is responsible for hysteresis in 
the stress-strain behaviour of M06S3I6. This interpreta- 
tion is slightly different from Yang et ali£ where the initial 
atomic configurations were subjected to random distor- 
tions, which of course mimic the entropy which helps to 
overcome the potential barriers in the metastable region. 

Although the Mo-S-Mo bonds in M06S3L3 nanowires 
are very deformable, the Young modulus of the un- 
strained energy minimum is only slightly smaller than 
Y of the second minimum with strained Mo-S-Mo bonds. 
The lack of these very deformable Mo-S-Mo bridges in 
MogSe is also the reason why the latter chains have a 
higher Young modulus. Surprisingly, the chains break 
in the metallic-covalent bonds between Mo in octahedra 
and not, as one would expect, in the polar covalent Mo- 
S bonds across the S bridges. This is the reason why 
all Mo nanowires, dressed with I or S anions and with 
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or without bridging atoms, have similar calculated ten- 
sile strengths. However, the much larger experimental 
Young's moduli reported by Kis et remain a puzzle 
at present. The M06S3I6 nanowires behave as quasi- one- 
dimensional conductors in the whole range of investigated 
strains. The conductivity involves the Mo octahedra and 
the bridging atoms and is extremely sensitive to strain, 
making this material very suitable for stain gauges. 

The fact that the experimentally reported shear mod- 
ulus appears to be an order of magnitude smaller than 
in SWCNTs suggests that ID effects will be even more 
pronounced than in SWCNT ropes, and the low-energy 
spectrum as E — > Ep is expected to be modified as a 
result. Very clean M0SI2; nanowires with good contacts 
may be expected to behave as ballistic quantum wires 
over lengths of several fim. On the other hand, with high- 
impedance contacts MoSI x nanowires may be thought of 
excellent candidates for the observation of Luttinger liq- 



uid behaviou»i2iiS. Of course, one-dimensional metallic 
systems such as M06S3I6 are subject to strong localisa- 
tion effects. Defects on a molecular wire can effectively 
stop and localise electrons, whereupon hopping between 
molecular chains becomes relevant. The pronounced ID 
nature of the nanowires makes them a uniquely versa- 
tile and user-friendly system for the investigation of ID 
physics. 
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